Electronic thermoelectric power factor and Metal - Insulator transition in FeSb 2 
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We show that synthesis-induced Metal -Insulator transition (MIT) for electronic transport along 
the orthorombic c axis of FeSb2 single crystals has greatly enhanced electrical conductivity while 
keeping the thermopower at a relatively high level. By this means, the thermoelectric power factor is 
enhanced to a new record high S 2 <r ~ 8000 /iWK -2 cm -1 at 28 K. We find that the large thermopower 
in FeSb2 can be rationalized within the correlated electron model with two bands having large 
quasiparaticle disparity, whereas MIT is induced by subtle structural differences. The results in this 
work testify that correlated electrons can produce extreme power factor values. 

PACS numbers: 72. 20. Pa, 71.27.+a, 71.30.+h 



I. INTRODUCTION 

The efficiency of thermoelectric (TE) material at tem- 
perature T is evaluated by the figure of merit ZT = 
(S 2 /pn)T, where S is the thermopower, p is the electri- 
cal resistivity and k is the thermal conductivity. There 
are two distinct approaches to increasing the ZT: cither 
by k reduction or by power factor {S 2 / p) enhancement. 
Techniques such as alloying, " phonon glass electron crys- 
tal" (PGEC) approach, and nanostructure engineering 
have been used to reduce the phonon mean free path, re- 
duce lattice k and produce high Z, at or above the room 
temperature^^ Reduction of the phonon mean free path 
has limitations since it cannot be reduced below the inter- 
atomic spacing and therefore other mechanisms for high 
ZT are sought afterP Tuning the carrier density by dop- 
ing is often inadequate since lower p comes with higher 
carrier concentration that favors lower 

The most favorable TE electronic structure is the one 
that has a resonance in the density of states centered 
about 2-3 kgT away from the Fermi energy (e F )P Kondo 
Insulators (KI) represent a close approximation of such 
an ideal case. In a strongly correlated electron system, 
and, in particular, in a KI, localized / or d states hy- 
bridize with conduction electron states leading to the 
formation of a small hybridization gap. The density of 
states just below and just above the hybridization gap 
becomes very large. The thermopower is very sensitive 
to variations in density of states in the vicinity of the 
e F , hence very large absolute values of |5| > 100 /iV/K 
can be expected and are reported in KlP The materials 
include not only rare-earth based compounds but also 
FeSiPKU 

FeSb2 crystallizes in Pnnm orthorhombic structure 
and has been characterized as an example of strongly 
correlated non-cubic Kondo insulator-like material with 
3d ionsP^HII] Similar to FeSi, heavy fermion states were 
discovered in FeSb2 by doping-induced metallizationP^Sl 
Colossal values of thermopower up to ~ 45 mV/K at 10 
K and a record high thermoelectric power factor (TPF) 
of ~ 2300 /iW/K 2 cm were observed This is two orders 



of magnitude larger than in best Bi2Te3-based materi- 
als and one order of magnitude larger than any reported 
value in correlated electron systems. Crystals in Refs. 19 
and 21 exhibit semiconductor behavior and small resis- 
tivity anisotropy with relatively high resistivity for cur- 
rent directed along all three principal crystalline axes. In 
contrast, single crystals used in Kondo insulator studies 
showed high anisotropy of the resistivity: electric trans- 
port along the a and b axes is semiconducting, while 
the c axis is metallic at high temperatur es and e xhibits 
metal-insulator transition (MIT) at 40 K) 16 ' 17 ' 22 ' 

Here we investigated the TPF in two FeSb2 crystals, 
with (crystal 1) and without (crystal 2) MIT. We show 
that large S enhancement can be attributed to electronic 
correlations in multiple charge and heat carrying bands. 
We also provide evidence for structural origin of MIT in 
correlated electron bands that reduces p by several orders 
of magnitude around 30 K in crystal 1 and results in a 
new record high TPF. 



II. EXPERIMENT 

Single crystals of FeSb2 with MIT have been prepared 
as in Refs. 16 and 17 by decanting at 650°C. Single crys- 
tals of FeSb2 without MIT have been prepared by the 
method described in Refs. 19 and 21, i.e., by decanting 
at 690° C, after a cool down from high temperatures to 
640° C. Crystals were oriented using Laue camera and cut 
into 0.6x0.5x4.5 (axbxc) mm 3 samples for two probe 
thermopower measurements with two ends soldered on 
disk-shape leads using indium. Resistivity p(T) was de- 
termined by a standard four-point ac method. Heat and 
current transport along the orthorombic c axis were mea- 
sured using the Quantum Design PPMS platform. A 
magnetic field was applied along the crystallographic a 
axis. Sample dimensions were measured with an optical 
microscope Nikon SMZ-800 with 10 /im resolution. Con- 
sequently the relative errors on the electrical resistivity 
and thermal conductivity are 4% and for the Hall coeffi- 
cient 2%. Since the Seebeck coefficient does not depend 
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FIG. 1. (a) Comparison of c-axis resistivity for two crystals. 
In 9 T p(T) in crystal 1 above 80 K is consistent with semicon- 
ducting thermally-activated transport with an activation gap 
of 120 K. (b),(c) Temperature regime(s) where p(T) is dom- 
inated by thermally activated transport. The arrow in (b) 
marks the MIT in crystal 1, where d[ln(p)]/dT changes sign. 
Broken lines in both panels are linear fits to activation-type 
behavior. 

on the sample geometry the main source of error is the 
sample uniformity and the accuracy of crystal orientation 
which introduces the measurement error of up to 5%. 

The atomic pair distribution (PDF) method, based 
on the total scattering approach, yields structural in- 
formation on different length scalesP^l X-ray scattering 
experiments for PDF analysis were carried out at the 
11-ID-C beamline of the Advanced Photon Source us- 
ing high energy beam (E = 114.82 keV, A = 0.108 A, 
0.5 mm x 0.5 mm size). Both PDF, G(r), and its alge- 
braic relative, radial distribution function (RDF), R(r), 
were considered in this study. Experimental PDF, G(r) 
is obtained from the measured reduced total scatter- 
ing structure factor, F(Q)=Q[S(Q)-1], via sine Fourier 

oo 

transform G(r) = G{r) = (2/tt) J {F(Q)Sin(Qr)dQ). 

o 

In practice, the upper limit of integration is some finite 
value Q m ax- RDF, R(r) is obtained from G(r) through: 
R(r)—rG(r) + 4irr 2 po where po is the average number 
density. Experimental setup for total scattering x-ray ex- 
periments utilized a Cryo Industries of America cryostat 
and Perkin-Elmer image plate detector. Finely ground 
samples in cylindrical polyimide capillaries were placed in 
a low-temperature sample changer, and the data for the 
two samples were successively collected for 4 min. at each 
temperature in probed range between 5 and 300 K. PDFs 
were obtained up to Q ma x = 26 A -1 momentum transfer 
using standard protocols, and intermediate length scale 
structure modeled over (1.7-45.0) A" 1 range with Pnnm 
structural model using the program PDFGUlpB Prelimi- 
nary reference neutron total scattering based PDF's were 



obtained using a time-of-flight HIPD instrument at Los 
Alamos Neutron Scatttering Center. 

III. RESULTS AND DISCUSSIONS 

The p c for crystal 1 is metallic down to th e onset of 
MIT (Tmit), as opposed to crystal 2 [Fig. l(a)]! 1B ' 17 ' 21 ' ^ 
Above 80 K and up to 300 K, p c (T) of crystal 2 is insen- 
sitive to magnetic field and is semiconducting. Interest- 
ingly, the logarithmic derivative, d[lnp c (T)]/dT reveals 
an underdeveloped anomaly at Tmit in crystal 2 which 
points to the intrinsic nature of this temperature scale 
[Fig. l(b,c)]. Arrhenius analysis assuming thermally ac- 
tivated behavior leads to a number of distinct energy 
scales valid for limited subranges of this temperature in- 
terval (Table I), allowing only approximate estimation 
of the gap values. The absolute value \d[ln(p(T))]/dT\ in 
the metallic and insulating states is very similar, which is 
reminiscent of the universality observed near the MIT in 
low-dimensional electron gas systems.^ Below 5 K elec- 
trical transport is dominated by extrinsic impurity states 
and is governed by a very small activation gap, similar 
for both crystals. Distinct energy scales and indirect gaps 
are in agreement wit h optical spectroscopy studies and 
ab initio calculations! 18 * 21 * 26 * 27 * The large activation gap 
observed above 80 K in crystal 2 is not observed in the 
metallic phase of crystal 1, but its size argues against 
contribution of impurity states. 

The Hall resistivity p xy (B) (current flowing along the c 
axis) for both crystals [Fig 2 (a,b)] is nonlinear, confirm- 
ing the presence of two carrier bands. In the absence of 
skew scattering we proceed to the analysis of the data in 
the two-band picture. In a two-carrier system^ the Hall 
coefficient is R H = p xy /H = p (a 2 + foH 2 )/(l + (3 3 H 2 ) 
where a 2 = /i^i+/2^2, /?2 = (/iM2+/2Mi)a*iA*2 and/3 3 = 
(A>"2 + /2A*i) 2 where p = p(B = 0), / 4 = \niHi\fE\niPi\ 
is the f factor, rii and pi are individual carrier band 
concentrations and mobilities. The agreement with the 
model is excellent and we obtain carrier band concen- 
trations and mobilities for crystal 2 while also showing 
data for crystal 1 for comparison.^ l n the model " car- 



TABLE I. Approximate (see text) energy gap from elec- 
tronic transport, the ratio of quasiparticle weights and scat- 
tering times in the valence and conduction band ZY = 
([(Z t , 2 r„ 5/2 )/r„]/[(Z c 2 r c 5/2 )/r c ])and e F ~\n{p c /p v ) for a crys- 
tal with MIT (top) and a crystal with no MIT (bottom 
row). The huge difference of fitting parameters between crys- 
tals 1 and 2 may be rationalized by the fact that the mobility 
difference of two bands for crystals with MIT is 10 8 -10 9 while 
the one for crystals without MIT is ~10, whereas both show 
similar large Seebeck coefficients. 
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FIG. 2. The p xy (a) and Hall constant (b) for crystal 2. Hall 
mobility (c) where the hole (electron) carriers are denoted 
by + (-) an d carrier concentrations of high and low mobility 
carriers (d) as a function of temperature for both crystals. 



rier" denotes a set of carriers with identical mobility as- 
sociated with only one energy and/or one degenerate en- 
ergy level and is different from the conventional electron 
or hole carriers that correspond to a continuous energy 
band. The low mobility carriers in crystal 2 are hole type 
(/x(2) l ~ 10 cm 2 /Vs) and they are an order of magnitude 
less mobile when compared to hole carriers in crystal 1 
at 300 K. Note that hole carriers in crystal 1 at room 
temperature constitute high mobility band (p(1)h ~ 10 2 
cm 2 /Vs). Yet, just like /J,(1)h, the h(2)l band exhibits 
sign change below 50 K, having identical temperature 
dependence when compared to with up to three 

orders of magnitude lower nominal values. The /x(2)jj 
carriers are electron type in crystal 2. They show little 
change in magnitude and do not change sign. The most 
striking difference between crystals 1 and 2 is the ratio 
of nominal mobility values. The large mobility difference 
is absent in crystal 2. Moreover, there is about or less 
than one order of magnitude between /j.(2)h and /i(2)^ 
in the high thermopower region from 10 to 40 K. The 
low and high carrier concentrations in both crystals have 
nearly identical values, suggesting a compensated nature 
of electronic transport at all temperatures. The carrier 
concentrations ri2(H) and n%(L) are rather close at 300 
K but they are both about an order of magnitude lower 
and have the same temperature dependence as carrier 



concentrations ni(H) and n\{L) in crystal 1. Carrier 
concentrations in crystals 1 and 2 are nearly identical 
in the temperature region of high thermopower, whereas 
in the impurity regime below 10 K crystal 2 (no MIT) 
shows much higher carrier concentrations. Values of car- 
rier mobilities for crystal 2 (^10 -3 cm 2 /Vs) are in gen- 
eral agreement with mobility values at low temperatures 
obtained using a two-band model (4 - 25 K) on crystals 
with no MITP 

Using carrier concentrations and Hall mobilities from 
a large number of measured p xy {B) isotherms we ob- 
tain S(T) for both crystals in the framework of a two- 
band non-interacting semiconductop 3 ^" 3 ^ model where 
thermopower S — {S e a e + Sh(Th)/((Je + °7i) and 



S e ,h = (k B /e)[ 



where Fj(£) 
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The scattering exponent 



S l+exp(x — £) ' 

s represents the energy dependence of the relaxation 
time t=to£ s , £=£F/(k,BT) are the reduced Fermi en- 
ergies for electrons and holes £ e and ^ (as measured 
from the bottom of the conduction band for electrons 
and from the top of the valence band for holes) and 
= h 2 /(2m*)(N/V)( 2 / 3 )(3/87r)( 2 / 3 ). Hence, for semimet- 
als e Fe = - {(Fh + eo) and £ e = - £/, - e /(k B T) where 
eo is the overlap energy, and for semiconductors ep e = 
- f-Fh + e 3 and £ e = - £ h - e /(k B T) where e g is the 
energy gap for semiconductors. In this model we used 
value s = -1/2 assuming that acoustic phonon scat- 
tering is dominant at all temperatures and m* = m e . 
The two-band non-interacting semiconductor model ad- 
equately describes crystal 2 above 50 K and below 5 K, 
whereas it fails to explain |S(T)| for crystal 1 in its metal- 
lic state (Fig. 3). While the similar shape of calculated 
and measured thermopower for both crystals between 10 
and 40 K in the noninteracting model argues in favor 
of the two-band approach, such calculation does not ex- 
plain the amplification of |5(T)| observed in the (10-30) 
K range. The prime suspect for discrepancy are strong 
correlations which can have a significant impact on car- 
riers in 3d bands. 

In a correlated electronic system with two bands sep- 
arated by a gap, S depends not only on the gap size but 
also on the anisotropy or asymmetry in the transport 
functionP 
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where tp is the chemical potential, A is the gap 
and SX=(X C -X V ) /(X c +X v ) is the asymetry parame- 
ter. The asymmetry parameter carries the infor- 
mation on band specific correlation strengths since 



X r ,.v — 



*(6/2).-fc, 



e P^/T^vTKq where Z c _ v are quasipar- 



ticle weights, r CjlI are scattering amplitudes and m* „ 
is the effective mass of carriers in conduction and va- 
lence bands. In addition, in real materials the \S(T)\ max 
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rapidly diminishes with the increase in impurity con- 
centration due to the increase of scattering amplitude, 
as recently observed in FeSb^f^^ Bandwidth narrow- 
ings to* v influence S(T) not only through the asymme- 
try parameter S\ but also via chemical potential since 
e F =(3k s T/4)ln(m*/m*)P Finally, the presence of im- 
purities can have considerable effect on both \S(T)\ max 
and S(T). The S(T) magnitude at a given tempera- 
ture depends on the ability of impurity carriers to put 
the Fermi level in the optimal position for thermopower 
enhancement P 

Both crystals have similar S(T) and nearly identical 
magnetothermopower MT = [S(9T)-S(0)]/S(0)J in the 
region of high S (Fig. 3). In contrast to differences 
in p(T), thermopower S(T) for heat flow along the c- 
axis is rather similar above 120 K and below 8 K. Ther- 
mopower changes sign from positive to negative above 
100 K in both crystals, indicating the presence of two 
carrier types. Since to* v =er CjV / /j, C;V where t c>v is the 
scattering time and /i c . v is the mobility in conduction 
and valence bands, we use mobility values for individual 
bands obtained from the fits of the Rjj in the two-band 
model. The value of Z 2 t 5 ^ 2 /T then becomes a fit param- 
eter for each carrier band, in addition to the chemical 
potential ep . Fits to the above equation for a fixed value 
of the gap which corresponds to the correlated electron 
temperature region (10-30) K are shown in Fig. 3(a) 
as red (crystal 1) and blue (crystal 2) solid lines. The 
best fits are obtained for the gap values of ~ (15-20) 
meV (Table I), suggesting that the large enhancement 
of S is due to the strong electro nic co rrelations associ- 
ated with a smaller indirect gap! 18 ! 26 ! Note that this is 
within the extremal limits in the assymetric case when 
the chemical potential is near the edge of the valence 
band (|S(T)e|< A/T+5/2k B )P For large T or for simi- 
lar A c and A„ (symmetric multiband effects) S —>0. This 
is observed: S sign change is at T = (120± 4) K, for both 
crystals 1 and 2. This implies that correlation effect on 
S(T) vane at high temperatures. 

The phonon drag is unlikely to have a significant 
contribution to S since isostructural RuSb2 and FeAs2 
have larger re and much smaller values of \S(T)\ max due 
to differe nt tem perature dependence of S and Nernst 
coefficient PlIlMl Tne | S ( T )| rises below 40 K, reach- 
ing values of 1 mV/K (crystal 1) and 1.9 mV/K (crystal 
2) in the (10-20) K range (Fig. 4). The S 2 a of crystal 1 
(7800 jzWK _2 cm _1 ) is maximized at T max = 28 K, which 
is higher than in crystal 2 (T max = 20 K) and close to 
Tm/t- The maximum value we find is three times larger 
than the TPF previously observed^ in FeSb2 and occurs 
at 16 K higher temperature. 

The central finding of the above analysis is that large 
thermopower in both FeSb2 crystals can be rationalized 
within the correlated electron model with two bands hav- 
ing the large disparity of quasiparticle properties. This 
is best illustrated by the ratio of carrier mobilities in the 
conduction and valence bands. While the mobility of the 
valence band carriers in crystal 1 greatly exceeds that of 
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FIG. 3. (a) \S(T)\ for FeSb 2 crystals 1 (red) and 2 (blue 
squares). Red and blue balls show fits to the two band non- 
interacting semiconductor model. Fits to the correlated elec- 
tron model for crystal 1 are shown by red, purple and brown 
(crystal 1) and dark blue, light blue and violet (crystal 2) solid 
lines that correspond to increasing values of fixed energy gap 
A (see text). The ratio of individual band mobilities in a 
two-carrier model is shown by the red/orange (crystal 1) and 
dark/light blue squares (crystal 2). 




T(K) 



FIG. 4. (The low resistivity around MIT leads to a record 
high TPF. Crystal 1 has two orders of magnitude higher TPF 
between 8 and 100 K. 



the conduction carriers, they are practically immobilized 
in crystal 2 where conduction carriers have higher mobil- 
ity. A key observation here is that despite great disparity 
in absolute values, mobilities of the valence band carriers 
in two crystals show strikingly similar temperature de- 
pendence with a pronunced anomaly at MIT. Most im- 
portantly, in both cases there is a change in the nature of 
valence charge carriers, which are holes at temperatures 
above MIT and electrons at lower temperatures. As we 
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discuss later, this observation reveals the likely nature of 
the MIT and of its absence in crystal 2. 

Now we turn to the quasi-one-dimensional (quasi-lD) 
metallic conductivity^ and the MIT in crystal 1 which 
is the key for colossal TPF. For a large concentration 
of impurities the chemical potential can go into the va- 
lence/conduction band, producing metallic p{T)W^ While 
an order of magnitude disparity in the number of charge 
carriers observed in crystals 1 and 2 above ^100 K would 
be consistent with such a scenario, it cannot be recon- 
ciled with the closely compensated nature of the carrier 
content. Below 10 K, where conductivity is governed by 
impurities, there are more carriers in crystal 2 than in 
crystal 1. Finally, impurity bands cannot account for the 
quasi-lD metallic conductance. Band structure suggests 
that the likely origin of the quasi-lD transport is the non- 
bonding d xy band, where the overlaps of Fe d xy orbitals 
are along the chains of edge-sharing octahedra parallel 
to the c axis, with little or no overlaps between orbitals 
in different chains [Fig. 5 (a)]! 27 * 40 1 ^ Hall data shows 
that quasi- ID metallic conductance in sample 1 at tem- 
peratures above 40 K is provided by a small number of 
mobile holes, which implies a nearly filled valence band. 
Such a non-bonding band must be narrow and strongly 
correlated, and should be described by the ID Hubbard 
modelP^ The low-dimensional transport in such band is 
sensitive to disorder and MIT is expected at half filling. 
Can it be the origin of the MIT observed in crystal 1 that 
d xy band is depleted with the decreasing temperature, 
attaining half-filling at Tmjt? Could it then be that a 
small additional disorder in crystal 2 induces localization 
of the quasi- ID charge carriers? 

In order to answer this question, we performed the 
PDF analysis of the x-ray data taken on powder samples 
of crystals 1 and 2 (Figs. 5 and 6). The PDF provides 
insight into the crystal structure on short and interme- 
diate length-scales and allows quantifying subtle struc- 
tural features such as bond properties and local disor- 
der. Intrachain and interchain Fe-Sb-Fe angles, /?i(T) 
and ,82 (T) respectively, shown in Fig. 5(b), are very sim- 
ilar, with a small difference developing in f3\ above 100 
K. This may indicate charge redistribution involving d xy 
bands since the presence of Fe d xy charge will result in 
a somewhat larger fti angled Atomic displacement fac- 
tors (ADP) [Fig. 5(c)] of Fe atoms are rather anisotropic: 
mean-square atomic displacements along the a and the 
b axes, Uu and U22, are nearly equal for both crystals, 
while those along the c axis, ^33, are noticeably enhanced 
in crystal 2 at all temperatures, being 0.0084 A (crystal 
2) and 0.0054 A 2 (crystal 1) at 10 K. Enhanced ADPs 
are typical indication of disorder. In this case disorder is 
markedly ID, precisely along the c axis and significantly 
higher in crystal 2. 

It is clear from the difference curves (Fig. 6) in the 
x-ray RDF (Refs. 23 and 24) that the average local en- 
vironments on the length scale of about one Pnnm unit 
cell (reflected in the first several RDF peaks up to 4.5 A) 
are indistinguishable. On larger length scales structural 
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FIG. 5. (Structural unit of FeSt>2 (a) shows 2x2x2 Pnnm unit 
cell as seen from the top of the c-axis (left), and from the top 
of the a-axis (right). Shaded rectangles denote short Fe-Sb 
distances within the FeSb6 octahedra. Intrachain and inter- 
chain Fe-Sb-Fe angles are denoted as /3i and 02, respectively. 
Next-nearest-neighbour interchain Fe-Fe distance is indicated 
by the double arrow, (b) Temperature dependence of /3i and 
P2 for crystal 1 (red) and crystal 2 (blue), (c) ADP factors 
for Fe (anisotropic) and Sb (isotropic) for crystal 1 (red) and 
crystal 2 (blue). Vertical dashed lines in (b) and (c) at 40K 
and 100K indicate temperatures of insulator-metal transition 
and change of the nature of the charge carriers, respectively. 

features start to differ. The RDF peak at about 4.65 A, 
denoted by a vertical arrow in Fig. 6, is markedly broader 
and less intense for crystal 2 at all temperatures studied. 
This peak corresponds to the next-nearest-neighbor Fe- 
Fe distance [FeSbg interchain distance, double arrow in 
Figure 5(a)], with no other contributions. Such a discrep- 
ancy could have several origins arising in a difference in 
Fe-Fe bondlength distribution, amount of Fe in the struc- 
tures of the two crystals, charge state of Fe between the 
two crystals, or a combination of these effects. 

Low-dimensional p is very sensitive to disorder which, 
in the simplest picture, suppresses metallic state by in- 
ducing strong localization. More generally, the disorder 
can impact the d xy overlap and the band structure, the 
orbital character of the electronic states responsible for 
conduction, the occupancy of Fe d xy orbitals, and even 
the orbital-dependent Hubbard U interaction strength 
in a d xy quasi-lD band of itinerant states that forms 
with sufficient d xy overlap. As opposed to this "self- 
generated" impurity level arising from Fe d orbitals in 
crystals with MIT, extrinsic impurity levels cannot pro- 
duce metallic resistivity near or close to room tempera- 
ture as in crystal s with MIT and are significant only at 
low temperaturespSHSl 

Crystals with MIT were cooled slowly from 1000 °C to 
650 °C and were decanted at that temperature from the 
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FIG. 6. Experimental neutron RDF data (inset) at 20 K 
reveals no detectable difference, suggesting that the disorder 
observed in x-ray experiment probably originates from the 
charge sector. 

liquid Sb flux. On the other hand, crystals with no MIT 
were cooled to 640° C, closer to Sb solidification temper- 
ature and were subsequently decanted at 690° C. There 
are two main factors that could contribute to increased 
crystallographic disorder in crystals with no MIT. FeSb 2 
melts and decomposes at 738° C to FeSb and Sb. Hence, 
crystallization occurs over a rather narrow temperature 
window. Crystals decanted at 690° C are much closer to 
the melting point when compared to crystals decanted 
at 650° C. This could contribute to increased disorder in 
Fe-Sb chemical bonds. In addition, crystals cooled to 
640 °C are closer to solidification line and are likely to 
experience more stress from the flux. This is consistent 
with observed structural differences, suggesting that MIT 
in FeSb2 crystals is governed by subtle structural differ- 



ences tunable by synthesis procedure. Greatly increased 
conductivity near MIT, combined with significant elec- 
tronic TEP, leads to a new record high thermoelectric 
power factor. 



IV. CONCLUSION 

In summary, we report the highest known TPF in- 
duced by the synthesis-controlled MIT in the correlated 
electron semiconductor FeSb2- The large thermopower 
enhancement can be understood within the electronic 
model, whereas MIT likely originates in quasi-lD and 
strongly correlated narrow band of itinerant states very 
sensitive to disorder. This is further supported by the re- 
cent observation of quasi-lD magnetism in isostructural 
CrSb 2 P 
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